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Abstract 

We analytically prove the necessary and sufficient criterion 
for the full separability of three-qubit Greenberger-Horne- 
Zeilinger (GHZ) diagonal states. The corresponding en- 
tanglement is exactly calculable for some GHZ diagonal 
states and is tractable for the others using the relative 
entropy of entanglement. We show that the biseparable 
criterion and the genuine entanglement are determined 
only by the biggest GHZ diagonal element regardless of 
all the other smaller diagonal elements. We have com- 
pletely solved the entanglement problems of three-qubit 
GHZ diagonal states. 

PACS number(s): 03.67.Mn; 03.65.Ud 

Introduction. — Quantum entanglement, being a special 
form of quantum superposition, possesses structures and 
properties intrinsically different from any classical system. 
This difference depends on the number of particles or 
states involved in the entanglement. For example, the test 
of Bell inequality using two-qubit entangled states gives 
statistical prediction, while the three-qubit Greenberger- 
Horne-Zeilinger (GHZ) entangled states lead to a conffict 
with local realism non-statistically [T,T . Furthermore, the 
entangled states of more than two qubits are more com- 
plicated because of a complex structure due to different 
ways the qubits can be entangled [3^ . Therefore new char- 
acteristic methods for multi-qubit states are necessary to 
fully understand and interpret their quantum behaviors in 
quantum mechanics. 

Multipartite entanglement right now is at the core of 
quantum information, and provides a critical resource for 
quantum secret sharing, quantum error-correcting codes 
and quantum computation [H [S]. A potential route 
to quantum computer is the manipulation of electrons 
trapped in quantum-dot pairs [5]. Recently, the first step 
revealing the true scalability of spin-based quantum com- 
puting was taken by coherently manipulating three indi- 
vidual electron spins confined in neighboring quantum dots 
[7], where the exchange interactions between the spins 
are precisely controlled and entangled three-spin states 



are generated. To date, multipartite entanglement has 
been observed in ion traps [8] , photon polarization [9] , su- 
perconducting phase and circuit qubit systems |10) . and 
nitrogen- vacancy centers in diamond [11]. Due to detri- 
mental decoherence effects and imperfections in prepara- 
tion, the multipartite entangled states prepared are usually 
mixed, typically, so called GHZ diagonal states. With the 
experimental realizations, a general separable criterion for 
GHZ mixed states is important and desired. 

One of the key issues is to determine whether the pre- 
pared states are genuinely entangled or not entangled at 
all. Theoretical research has concentrated on the sepa- 
rability and biseparability that characterize the entangle- 
ment of the GHZ diagonal states [12l[l3[Tl[T5l[li[T7l[18]. 
The criterion for biseparability of three and four-qubit 
GHZ diagonal states have recently been proved in the form 
of an inequality involving several density matrix elements 
in the computational basis [13] . We will show that this cri- 
terion can be simplified to the largest GHZ diagonal com- 
ponent being equal or less than 1/2. This is a significant 
step that leads to the quantification of genuine entangle- 
ment. For the full separability of three-qubit GHZ diago- 
nal states, the sufficient criterion is proposed by directly 
constructing the fully separable states TT and the neces- 
sary criterion is proposed by the method of convex combi- 
nations [18]. Numerical results strongly indicate that the 
two criteria should coincide. However, it is inconclusive in 
the absence of an analytical proof. 

In this letter, we will prove the unification of the criteria 
for full separability. The full separability and biseparabil- 
ity provide a practical and accurate method for complete 
classification of three-qubit GHZ diagonal states, which 
now can be made into three categories: separable, entan- 
gled, and genuinely entangled states discriminated by full 
separability and biseparability, respectively. In addition to 
the complete entanglement structure of three-qubit GHZ 
diagonal states, we also find a new group of entangled 
states due to the cooperation of the off-diagonal elements 
in computational basis. This cooperative entanglement 
may not be detected by the commonly used positive partial 
transpose (PPT) criterion [19]. We will give the formula 



of genuine entanglement and reduce the quantification of 
entanglement to algebraic calculations. 

Necessary and sufficient condition of fully separable 
three- qubit GHZ states. — The three-qubit GHZ diagonal 
states take the form 

8 

p^Y.Pk\GHZk){GHZk\, (1) 
fe=i 

where the pk form a probability distribution. The 
GHZ state basis consists of eight vectors \GHZk) — 
-^{\Qx2X3) ± 11x2X3)), with Xi.Xi e {0,1} and Xi ^ xi. 
In the binary notation, fc — 1 = 0x2X3 for the '+' states 
and fc — 1 = 1x2X3 for the '-' states. Using Pauli matri- 
ces X, y, Z and 2x2 identity matrix /, the GHZ diagonal 
states can be written as 

8 

+KYYX + \^YXY + \^XYY\, (2) 

where tensor product symbols are omitted. 

The sufhcient condition of full separability of p is [18^ 



1 - |A_| - M > 0, 



(3) 



or 



l-IA-l-IAsl-IAel-IAyl-IAsI >0 (4) 

where A_ = min{A2 + A3 + A4, A2 - A3 - A4, -A2 + A3 - A4, 
— A2 — A3 + A4, }, and 

^ _ \J (A5A6 + A7A8)(A5A7 + A6A8)(A5A8 + AsAt) 

VAsAgAyAg 

The sufficient condition comes from an explicit construc- 
tion of (0 in a fully separable manner [T7] [T5] . 

Let us consider the 8x8 density matrix p with entries 
Pij in the basis |000) , |001) , . . . , |111) which are ordered 
in the canonical way. For GHZ diagonal states, the only 
nonzero elements of p are pa and Pi,9_i [i — 1, . . . , 8), and 
further we have pa = pg-i^g^i. 

The necessary condition for the full separability of GHZ 
diagonal state p [18] can be written as 



£(p,^) <C(^)k 



(6) 



where k = minjpjj (1 < i < 4)},]? = ((5, a,/?, 7) is a real 
vector, C{p, ]?) = Spis + ap27 + I3p3e + 7P54, and 



C(^) = sup[(5cos(a + 6 + c) + acos(a) + /3cos(5) 



-7 cos(c)]. 



(7) 



Here a, 6, and c are the angles. The relationship among the 
density matrix entries and the parameters A^ of a GHZ 
diagonal state is a simple linear transformation; 

(A5, -Ae, -A7, -As) = 4(pi8,P36,P27,P54)ff2, (8) 

K = (l-|A_|)/8, (9) 



where Ho is the 4x4 Hadamard matrix. Consider the triv- 
ial case that ^ is a positive vector, dt) achieves its max- 
imum S + a + P + j when a = 6 = c = 0. If the off-diagonal 
elements of p are all positive, the left-hand side of inequal- 
ity ([6]) is the probability mixture of the off-diagonal ele- 
ments. If two of the off-diagonal elements pi8, P27, Pae, P54 
are negative, the state can be transformed to a state with 
all positive off-diagonal elements by local operations and 
classical communication (LOCC). For these cases, the in- 
equality ^ can be rewritten as 

max{|p,,9_,| , (1 < i < 4)} < min{p,„ (1 < i < 4)}. (10) 

It is just the PPT criterion. In the case of even number 
of negative pis, P27, Pse, P54, namely JlLs < 0, PPT 
criterion is necessary and sufficient for the full separability 

m- 

What left is to show that if the inequality (|3]) can be de- 
rived from (l6]) or vice versa when nLs K > 0. When one 
of the components of is negative (for definiteness, one 
may choose 7 < 0), the solutions have been obtained [18] 
for the angles a, 6, c to optimize C(^) for any given coef- 
ficients S,a,(3,j. The solutions fulfill the following equa- 
tions derived from ([7]) 



6 sin d 



-asm a 



-/3sin5 = I7I sine. 



(11) 



where d= a + b + c. Without loss of generality, we consider 
the case where pi,Si P3,6i P2,7 are non-negative (this is al- 
ways possible by LOGC), denote {x,y,z) — {S,a,l3)/ I7I , 
and define 



f{x,y,z) 



C{pM 



{xpis + yp36 



ZP27 - P5i)/K 



X COS d + y cos a + z cos b — cos c 



(12) 



For the outmost surface of the full separability set deter- 
mined by the necessary condition, we have the following 
equations 



dfix,y,z) 
dx 
f{x,y,z) 



df{x,y,z) df{x,y,z) 



dy 



dz 



= 0, (13) 
(14) 



The solution to equations (IT3l) is 



P18 = n cos d, P36 — K cos a, P27 ~ n cos &, ^54 = cos c, 

(15) 

where we have used equations ([TT]) , ([T4]) and the fact that 
d = a + b + c and a, 6, c, d are functions of x, y, z in equation 

(US- 

Substituting the solution (llSp into equation ([8]) to obtain 
Xi {i — 5, ...,8), and using of ([5]) and (jH]), we have 



p = l-|A_|, 



(16) 
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Figure 1: The border curve of full separability 



Hence the outmost surface of fully separable state set de- 
termined by the necessary condition takes the form of Eq. 
(IT51) . which is also the outmost surface of fully separable 
state set determine by the sufficient condition ([3]). In ei- 
ther odd or even number of negative parameters S, a, /3, 7, 
it is shown that the necessary criterion of full separabil- 
ity is also sufficient. Thus we obtain the necessary and 
sufficient criterion for fully separable of three-qubit GHZ 
diagonal states. 

Consider the case of pis = P36 = P27, we have the simple 
solution of 6 = a, c = —3a, d = —a [a € [0, 7r/3]) for the 
border curve of the full separability. More explicitly, we 
have the equation of border curve 



K K 



,Pl8 



(17) 



The border curve is shown in Fig.l. 

Entanglement: — Entanglement of a given entangled 
state g can be quantified by its relative entropy of en- 
tanglement (REE) as, E — muiaeSEP S{g\\a), here SEP 
is the fully separable state set. The relative entropy is 
defined as, ^(gllcr) =: Trg\og2 g — Trp logjcr. The sepa- 
rable state achieves the minimal relative entropy for state 
is called its closest state. Follow the same reason given 
in ref f2U|, the closest state for a three-qubit GHZ diag- 
onal state p must also be a GHZ diagonal state in the 
form oi <T — J2k=i Ik \GHZk) {GHZk\, where the qk form 
a probability distribution. As the result, the REE can be 
calculated as, 

£; = Vp,iog2(^). 



fc=i 



qk 



(18) 



The closest state a is at the boundary of the fully sepa- 
rable states due to the following reasoning: Suppose a' is 
an inner state within the fully separable state set. The 
entangled state p is at the outside. A link between a' and 



p should intersect with the boundary of the fully separa- 
ble state set. We have a border state a = Act' + (1 — \)p 
for some A S (0,1). Then S'(pllcr) < XS{p\\a' ) + (1 - A) 
S{p\\p) = XS{p\\a') < S{p\\a'). The first inequality 
comes from the convexity of negative logarithmic func- 
tion. Hence an inner separable state is not a closest 
state. The parameters of a closest state can be written as 
qk = Sk + KcCos9k and qg-k = Sk-KcC0s9k for fc = 1, • • -4, 
where Sk >0,Kc = ininfc{sfc}, 6 = (d, a, b, c). The free vari- 
ables are a, b, c and three of Sk- We have six free variables 
and six equations derived from the extremal of (fT8|) . REE 
is tractable. In some special cases, analytical solutions can 
be obtained as we will show below. 

If an entangled GHZ diagonal state p has the sym- 
metry of pi = P2 = P3 and pe = pr = Ps, we will 
prove that the closest state should have the same sym- 
metry of qi — q2 = 93 and qe = (?7 ~ qs- Suppose that 
a nonsymmetric state ui = J2^^^^qk\GHZk) {GHZk\ is 
the closest state of p. By cycling (71,92,93 and qG,qr,qs, 
we obtain <T2 and (73 which also are closest states of p. 
Let a = ELi9fc \GHZk) {GHZk\ = i(ai + ^2 + ^3), 
with = 92 = 93 = i(9i + 92 + 93), 96 = 9? = 98 = 
■5(96 + 97 + 98), then a is fully separable since it is a 
probability mixture of fully separable states. Note that 
|(9i + 92 + 93) > sy9i9293- We have S{p \\a) < S{p \\ai ), 
that is a contradiction. The symmetric closest state 
cr is with an = (T22 = 0-33 and (ctis, 0-36, 0-27, 0-54) = 
Kc(cos0,cos6',cos0, cos36'), where 9 G [0, f], and Kc = 
minjcTii, (T44}. Let CTh = ^ + "'44 = ^ ~ where 
^ G ' '^^ obtain the two optimal equations with 

variables 9, The equations are not analytically solvable 
in general. However, if the closest state is with Kc = |, 
exact solution can be obtained. The derivative of (ITSt on 
9 leads to 



(- -4pii)i + pi8 -P54 = 0, (19) 



where t = 2 cos 6*. We then properly choose 9 from the 
three solutions of ([T5|) . Note that the derivative of ^ may 
not exist at ^ = 0, we have to calculate the left and right 
derivatives, which are 



dS{p\\a] 



=0+ 



ci(6')[cos36'(cos36'- 8p54) +8p44 - 1], 
02(6*) [cos 6'(8pi8 - cos 61) + 1 - 8pii], 



where we have used equation (jl9p to simplify the expres- 

8 J fa\ 8(4cos^6(-l) 

and 02(9) - ,i„.3e 



sions. ci{9) 



are positive 



3 sii?" w 

factors for 6* e [0, f ]. 

Analytical solutions exist for the states with pn = pa = 
I when 9 € (f,f]- The entanglement can be obtained 
exactly. We verify that cos 9 < Spis, cos 30 > 8/954 are true 
for all the entangled states. The left derivative is negative 
since cos 30 < for e (f , f ] and the right derivative is 



3 



positive. When P54 — — pis, the solution of (fT9| is 6* = ^. 
The REE is 



E = 



1 



5l8 , 1 
— l0g2 



-'18 



1 



1 + 75 



— l0g2 , _ j_ 
^ V2 



Analytical solutions also exist for states with pig = pn, 
P54 = -P44- Then pi = 2pii, = 2p44, ^4=^3 = 0. 
We have four kinds of candidate closest states: (i) Kc = \, 
(ii) Kc = (Ti l < |,(iii) Kc = cr4_4 < i, (iv) PPT boundary 
state. Here we consider > pi (The case ps < pi can be 
solved similarly). For the first closest state candidate, the 
relative entropy is simply 



S{p\\a) ^ 3pi log2 



Pi 



1 + cos ( 



P5 l0g2 



P5 



1 — COS 39 



(20) 



dS{p\\a) _ 



de 



= 0, we have (2 cos-^ 6* + A cose*- 1)(2 cos( 



From 

1) = 0, where A = ^aiiPi, The solution is cos 6* — 

' ^ P5+P1 

|;(V8 + A^ — A). To check that the closest state locates at 
Kc = |, we prove that the right derivative is always posi- 
tive and the left derivative is negative only when p\> p^ — 
(3 + cos 30) « 0.1718. Hence candidate (i) is the solution 
whenpi > po- For pi < po, we consider candidate (ii). We 
have S{p ||a) - 3pi loga ^-^^^^ +P5 logs i^.^o^^cosae) ' 



Then = 



gives - 2(3+ros3e) • We have 



5(p||a)-l + 3pi log2 



(3 + cos 36*) 
3(1 



cos ( 



(21) 



The optimal equation is 4 cos^ 6* + 6 cos^ — 'i = 0, the 
solution is cos 6* = i(\/2 + \/3+ \/2 - ^3 - 1) w 0.5979. 
The condition < ^ is equivalent to pi < po- The candi- 
date (iii) does not give rise to a further small value of the 
relative entropy because ^^jfj"^^ < for all G [0, f ]. The 
candidate (iv) can also be removed. 

Genuine entanglement: — The necessary and sufficient 
criterion has been proven for biseparability of a three-qubit 
GHZ diagonal state fTS'. It has the form of 



Also criterion (|23p and Eq. (|24p are true for four qubit 
GHZ state as the necessary and sufficient criterion has 
been proven [T3] . The mixture of N-qubit GHZ state and 
white noise is = {1 - p) \GHZN) {GHZN\ + 

pl/2^ . The state is genuinely entangled iff < P < 
1/[2(1 - 2-^)]. The probability of \GHZN) component 
is pi = (1 - p) + p/2^ = 1 - p(l - 2-^^) > i. The bisep- 
arable (yet inseparable under bipartitions) condition can 
also be written in the form of inequality (j23p . Therefore 
its genuine REE can be expressed as Ea. (p4)) . 

Summary: — For three-qubit GHZ diagonal state, the 
fully separable criterion has been strictly proven to be nec- 
essary and sufficient. We find the exact boundary states 
for the fully separable state set. The free variable solu- 
tion of boundary states make the calculation of the rela- 
tive entropy of entanglement easy. The relative entropy 
of entanglement is exactly calculated for the symmetric 
states p { pi = p2 = P3,pe = pj = ps) of type (i) pu = | 
{i — 1, ■ ■ ■ 8) (diagonal elements are equal in computational 
basis) and type (ii) p^ ~ ps = 0. The closest states a can 
be either with ct^^ = i (j = 1, • • • 8) or not for both type 
(i) and type (ii) entangled states. There are many sub- 
tleties in obtaining the closest states due to the possible 
non-existence of the derivative of the relative entropy. We 
also give the genuine entanglement of GHZ diagonal states 
in terms of the relative entropy of entanglement. The gen- 
uine entanglement is determined by the biggest GHZ diag- 
onal component only. The genuine entanglement formula 
obtained is easily extended to N particle GHZ diagonal 
states. The fully separable criterion can be applied as a 
necessary criterion for the separability of any three-qubit 
state by filtering it to a GHZ diagonal state. Further works 
on the separability and entanglement of more than three 
qubit GHZ diagonal states are desirable. 

XYC and LZ J thank the support of the National Natural 
Science Foundation of China (Grant No. 60972071). 
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The criterion can be rewritten as 



max{pj} < -, 



(22) 



(23) 



since for three-qubit GHZ diagonal state, we have pa = 
pg_j,g_j, and pi = -I- Ipi^sI • When pi > i, it is 
followed from the reasoning in the entanglement of the 
Bell diagonal state [20] that the state is not bisepara- 
ble. The genuine entanglement measured by REE is 

Egenmneip) = mi'O-aeBISEP S {p ||cr), whcrC BISEP IS the 

biseparable state set. Hence 



E„ 



,ip) = l+pilog2Pi + (l-pi)log2(l-pi). (24) 
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